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Abstract 

We introduce a two-dimensional (2D) discrete nonlinear Schrodinger (DNLS) equation with self- 
attractive cubic nonlinearity in a rotating reference frame. The model applies to a Bose-Einstein 
condensate stirred by a rotating strong optical lattice, or light propagation in a twisted bundle 
of nonlinear fibers. Two species of localized states are constructed: off-axis fundamental solitons 
(FSs), placed at distance R from the rotation pivot, and on-axis (R = 0) vortex solitons (VSs), 
with vorticities S = 1 and 2. At a fixed value of rotation frequency f2, a stability interval for the 
FSs is found in terms of the lattice coupling constant C, < C < C cr (R), with monotonically 

decreasing C CV (R). VSs with S = 1 have a stability interval, C cr (^) < C < C cr (^), which 

(s=i) 

exists for below a certain critical value, Ocr • This implies that the VSs with 5 = 1 are 
destabilized in the weak-coupling limit by the rotation. On the contrary, VSs with S = 2, that are 
known to be unstable in the standard DNLS equation, with = 0, are stabilized by the rotation in 

(5=2) (5=2) 

region < C < CS with CS 1 

growing as a function of Q.. Quadrupole and octupole on-axis 
solitons are considered too, their stability regions being weakly affected by ^ 0. 

PACS numbers: 03.75;42.65.Tg;05.45;42.70.Qs 
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I. INTRODUCTION 



Discrete dynamical systems represented by nonlinear lattices in one, two, and three di- 
mensions constitute a class of models which are of fundamental interest by themselves, and, 
simultaneously, they find applications of paramount importance in various fields of physics. 
One such example is known in nonlinear optics, where the one-dimensional (ID) discrete 
nonlinear Schrodinger (DNLS) equation was predicted by Christodoulides and Joseph to sup- 
port fundamental discrete solitons [l| with an even profile. Such nonlinear structures were 
later created experimentally in an array of semiconductor waveguides Subsequently, 
stable odd solitons, alias twisted localized modes, were predicted and studied in detail in 

n 

the same model [3[, as well as in an array of photorefractive waveguides with photovoltaic 
nonlinearity 0]. 

Another experimental realization of dynamical lattices in the optical domain is possible in 
photorefractive crystals, where a quasi-discrete setting can be induced by counterpropagating 
laser beams illuminating the crystal in the normal polarization, while the probe beam, which 
can sustain solitary waves, is launched in the extraordinary polarization. The difference 
from the array of waveguides fabricated in silica or in a semiconductor material is that the 
photorefractive nonlinearity is saturable, rather than cubic. This method of the creation of 
photonic lattices was proposed in Ref. [5[, and results obtained by means of the technique 
were reviewed in Refs. [6). In particular, both fundamental and twisted solitons in the ID 
lattice were reported in Ref. 7fl, and fundamental solitons (FSs) in the 2D lattice were 

n n 

created too [8[], as well as 2D vortex solitons (VSs) in the same setting [9]. 

Recently, the progress in the technology of writing permanent arrays of channels in silica 
slabs has made it possible to create 2D waveguiding lattices with the cubic nonlinearity, 
which emulate a bundle of nonlinear optical fibers with linear coupling in the transverse plane 

n 

[lyj • In particular, spatial lattice solitons of the surface and corner types were reported in this 



setting 111 (surface solitons were reported too in the 2D photonic lattice in a photorefractive 
crystal PL2j|)- Thus, ID and 2D discrete dynamical models have a potential for further 
applications to nonlinear optical media of various types. 

Another natural realization for the lattice systems is provided by a Bose-Einstein conden- 
sate (BEC) trapped in an optical lattice (OL). If the OL is strong enough, the underlying 
Gross-Pitaevskii equation (GPE) for the wave function in the continuum may be approxi- 
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mated by its DNLS counterpart [3, Q; the relevance of the discrete model in this setting 
has also been demonstrated experimentally in ID by testing its predictions experimentally, 

n n 

see, e.g., Ref. [15j. Higher-dimensional OLs can be easily created too [161] . Thus, the relevant 
DNLS equation may be one-, two-, or three-dimensional; in particular, various species of 3D 
discrete solitons, including those with intrinsic vorticity, have been predicted in this setting 



17j , and their stability has been systematically analyzed [18] . Still another implementation 



of the DNLS lattice in the space of any dimension, from one to three, is possible in terms of 
a crystal of microcavities trapping photons 19] or polaritons j^oj]. 

Discrete solitons of various kinds have been studied in detail theoretically in ID, 2D, and 
3D versions of the DNLS equation, see an earlier review 21] and the more recent works 
mentioned above. As mentioned above, some of these solitons have been created experimen- 
tally in optical media equipped with fabricated or photoinduced lattices. All these localized 
states have their counterparts in continuum models with periodic potentials that emulate 
the lattices. In particular, 2D solitons of both the fundamental and vortex types, which are 
unstable in uniform continua with the cubic self-focusing nonlinearity can be readily stabi- 
lized by the periodic OL potential 22J; for the stabilization of FSs, a quasi-lD potential is 
sufficient, instead of its full 2D counterpart 23] . Vortices are unstable too in the uniform 
space with the saturab.e nonlinearity Q, in which case they can a.so be stabilised by the 
periodic potential [9J. The relevance of these stabilization mechanisms in the continuum was 
also demonstrated for higher-order vortex solitons, and so-called supervortices, i.e., arrays 
of compact vortices with global vorticity imprinted onto the array, under both cubic and 
saturable nonlinearities [251 ] . 

Recently, it was shown that 2D solitons obeying the GPE in the 2D continuum can 



also be supported by a rotating OL [26j, 



271 ] . These solitons may be fully localized (spot- 



shaped) solutions to the equation with the self- focusing/ attractive cubic nonlinearity, placed 
at some distance from the rotation pivot and revolving in sync with the holding 2D lattice. 
In particular, the soliton can be placed at a local minimum of the rotating potential, while 
the pivot is set at a local maximum. These co-rotating strongly localized solitons are stable 
provided that the rotation frequency, Q, does not exceed a critical value, (^ cr ) mm - In the 
same model, but with a rapidly rotating OL, stable ring-shaped solitons (with zero vorticity), 
i.e., objects localized along the radius but delocalized in the azimuthal direction, have been 
found too, for Q exceeding another critical value, (^ C r) max - Note that the model does not 
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support any stable pattern in interval (^ cr ) min < < (^ cr ) max |26|. On the other hand, 
stable ring-shaped states, with both zero and nonzero vorticity, have been found in the 
model with the repulsive cubic nonlinearity and rotating quasi-one-dimensional (periodic) 
potential, if Q exceeds a respective critical value. Obviously, the latter model does not give 
rise to any localized state in the absence of the rotation. 

A rotating OL can be easily implemented in BEC experiments 



28 



then, if the OL is 



strong enough, it is natural to approximate the GPE in the co-rotating reference frame by 
an appropriate variety of the 2D DNLS equation. In addition to that, such a model may 
also describe the light propagation in a twisted bundle of nonlinear optical fibers, linearly 
coupled in the transverse plane by tunneling of light between the fiber cores. The objective 
of the present work is to introduce a model of the rotating discrete lattice, and find stable 
discrete solitons in it, both FSs and VSs. 

The paper is organized as follows. In Sec. II, we formulate the model, taking the un- 
derlying GPE in the reference frame co-rotating with the OL, and replacing the continuum 
equation by its discrete version corresponding to a strong periodic potential. Discrete FSs 
are considered in Sec. III. We construct the solutions starting from the anti-continuum 
limit, which corresponds to zero value of the coupling constant accounting for the linear 
interaction between neighboring sites of the discrete lattice, C = 0. A family of FS solutions 
is constructed by continuation in C; their stability is examined by computation of eigenfre- 
quencies for infinitesimal perturbations around the soliton, and verified by direct simulations 
of the evolution of perturbed FSs. It is found that the FS, with its center located at distance 
R from the rotation pivot, is stable within an interval < C < C CI (R), with C cr decaying as 
a function of R. Section III also includes a simple analytical approximation, which makes 
it possible to explain the decrease of C cr with the growth of R. In Sec. IV, we consider 
localized vortices (VSs), whose center coincides with the rotation pivot. For the VS with 
vorticity S = 1, the stability region is found to be Ccr < C < Ccr , provided that the 
rotation frequency, Q, is smaller than a critical value, f2 cr (the stability interval shrinks to 
nil at Q = fi cr ). Vortices with S = 2 are considered too. While in the ordinary (nonrotat- 
ing) DNLS model, with Q = 0, all VSs of the latter type are unstable 29|, we demonstrate 

(5=2) (S=2) 

that the rotation opens a stability window for them, < C < C cr , with C cr growing 
as a function of Q. Direct numerical simulations are also used to illustrate the dynamical 
evolution of FSs and VSs when they are unstable. Results obtained in this work and related 
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open problems are summarized in Sec. V. 



II. THE MODEL 

The starting point is the normalized 2D GPE, which includes the potential in the form 
of an OL rotating at angular velocity Q, and thus stirring a "pancake" -shaped (quasi-flat) 
Bose-Einstein condensate (BEC) trapped in a narrow gap between two strongly repelling 



optical sheets. Unlike the analysis performed in Refs. [26| and [27|], where simulations 
were run in the laboratory reference frame, here we write the GPE in the reference frame 
co-rotating with the lattice, hence the potential does not contain explicit time dependence: 

%d ~ti = ~ Q y2 + nLz ) ^ ~ e [c ° s ^ (x _ e)) + c ° s(A: (y ~ v))] ^ + (i) 

Here, L z = i(xd y — yd x ) = id$ is the operator of the z-component of the orbital momentum 
(8 is the polar angle), and a determines the sign of the interaction, attractive (a = — 1) or 
repulsive (a = +1). Constants £ and v determine a possible shift of the lattice with respect 
to the rotation pivot. Note that Eq. ([I]) does not contain any additional trapping potential, 
as we are interested in solutions localized under the action of the OL. 

It is more convenient to shift the origin of the Cartesian coordinates to a lattice node, 
thus replacing Eq. (CQ) by the translated form: 



.dip 



tp — e [cos (kx) + cos(ky)] ip + cr\ip\ 2 ip. (2) 



As shown in a general form in Ref. [14j . a discrete model, which corresponds to the limit of a 
very deep OL, can be derived from the underlying GPE in the tight-binding approximation. 
Eventually, it amounts to a straightforward discretization of the GPE. Thus, the discrete 
counterpart of Eq. ([2]) is 

i d ^' U = -^{{i } m+l,n+i>m-l,n+i>m,n+l+i>m,n-l — 4'lp m ,n) 

where (m,n) are discrete coordinates, and C > is the corresponding coupling constant, 
which accounts for the linear tunneling of atoms between BEC droplets trapped in deep 
nodes of the lattice. As mentioned above, Eq. may also describe a twisted bundle of 
nonlinear optical fibers linearly coupled by the tunneling of light in the transverse plane, 



(ra, n). In that case, t is the propagation distance along the fiber, and only a = — 1, i.e., the 
attractive/focusing nonlinearity) is the relevant choice. 
Equation conserves the norm and Hamiltonian, 

iV = ^|^ m ,n| 2 , (4) 

m,n 



f C 1 



1 4 



-— [{m + (VC,n (V»m,n+1 - ^m,n-l) ~ ^m,n fen+l ~ C*,n-l)) 
- (71 + V) (V4,n (V'm+l.n - ^m-l,n) ~ VV« (VC+l,n ~ ^m-l,J)] } • ( 5 ) 

In addition to C, the discrete model contains three irreducible parameters: Q, which takes 
values < Q < oo, and the coordinates of the pivot displacement, (£, i>), which take values 
< £, v < I; P ms the sign parameter, a = ±1. As in the usual 2D DNLS equation (with 
Q = 0), values a = ±1 in Eq. ([3]) may be transformed into each other by the staggering 
transformation, ip m . n — > (— l) m+n ip mtn , therefore we fix a = —1 (self-attraction). 

Our first objective is to find stationary localized solutions to Eq. ([3]) in the form of FSs 
(fundamental solitons) and VSs (vortex solitons). To this end, we substitute the standing 
wave ansatz, ipm,n — e lAt( firn,n, where —A is the normalized chemical potential, in terms of 
the underlying BEC model; then, the stationary lattice field 4> m ,n obeys the equation 

C 

^4>ra,n = ~ (0m+l ,n ~t~ (pm—l,n > (Pm,n+1 ~T~ (Prn,n—1 ^(Pm,n) 
C 

+i — [(m + f ) (4> m ,n+l ~ 4>m,n-l) ~ (n + v) ((f) m +l,n ~ <j>m-\,n)\ + \4> m ,n\ 2 4>m,n (6) 

Note that solutions for 4> m ,n are complex, unless Q = 0. The second objective will be to 
examine the stability of the discrete solitons, assuming small perturbations in the form of 
o~i/Jm,n ~ exp (iAt + iXt), the onset of instability indicated by the emergence of Im(A) 7^ 0. 
The evolution of unstable solitons will be examined by means of direct simulations of Eq. 

To parameterize the soliton families, we fix the scales by setting A = 1; obviously, with 
C > 0, only positive A may give rise to localized solutions), while C will be varied. 
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III. FUNDAMENTAL SOLITONS 



The rotation makes the discrete lattice inhomogeneous, hence properties of solitons 
strongly depend of the location of their centers. Without the rotation, Eq. ([3D amounts 
to the ordinary 2D DNLS equation, in which various species of discrete solitons and their 
stability have been studied in detail. In particular, FSs, which are represented by real so- 
lutions, are stable at C < C cr = 2 A = 2 |2l|. The onset of their instability is accounted 
for by a pair of eigenfrequencies of small perturbations with finite imaginary and zero real 
parts, i.e., the instability leads to the exponential growth of perturbations. Accordingly, nu- 
merical simulations of the instability development demonstrate spontaneous transformation 

n 

of unstable FSs into lattice breathers [30]. In this section, we first report numerical results 
obtained for the stability of FSs in the model with Q ^ 0, and then present an analytical 
estimate that may explain numerical findings. 



A. Numerical results 



Our analysis aimed to determine the stability border for the FSs, C cr , for each set of 
values of the discrete coordinates of the soliton's center, {m ,n }. Here we present results 
for angular velocity Q = 0.1 and zero pivot displacement, £ = v = [in Ref. 26(, the rotation 
pivot was fixed at a local maximum of the potential in Eq. (j2J), which corresponds to setting 
^ — v — 1/2 in Eqs. ([3]) and (jBJ), and the center of the soliton trapped in the lattice was 
placed at a local minimum closest to the pivot, which would mean {mo, no} = {0,0}]. This 
choice makes it possible to explore the existence and stability of FSs in a clear form, while 
larger values of Q give rise to a resonance with linear lattice modes, leading to Wannier-Stark 
ladders and hybrid solitons [33j and making the continuation in C and identification of C cr 
difficult. We carried out the calculations on the lattice of size 21 x 21, since already for this 
case, the lattice was for all the considered cases sufficiently wider than the very localized FS 
structures of interest. To avoid effects of the boundaries, the range of the soliton-center's 
coordinates was restricted to \w,q\ , |rto| < 8. 

The FS solutions were looked for starting at point C = 0, i.e., at the so-called anti- 
continuum limit 2lj. In this limit, the FS is seeded by using a nonzero value of the field at 
a single point, the center of the FS, 4>m,n°^ = 5m,m ^n,n ; obviously, this expression satisfies 
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Eq. Q) with A = 1 and C = 0. After the branch of the FS solutions had been found by the 
continuation in C, its stability was quantified through the computation of eigenfrequencies 
of small perturbations, using the equation linearized about the FS. 

Figure [T] displays a typical example of the thus found dependence of eigenfrequencies 
A of small perturbations on the lattice coupling constant. As said above, the instability 
corresponds to Im(A) 7^ 0, for the FS with its center set at point (m = 3,n = 2); for 
comparison, the dependence of the instability growth rate, i.e., |Im(A)|, on C is also shown 
for the FS in the usual DNLS model (Q = 0) in the top panel of the figure. It is seen that 
the instability sets in at C = C cr = 1.70, which is smaller than the critical value in the 
ordinary model, Cc?~°^ = 2. Typical examples of stable and unstable FSs belonging to the 
family presented in Fig. [T] are displayed in Figs. [2] and [3J 

Results obtained for the FSs placed at different positions are summarized in Fig. 2J in the 
form of dependences of C cr on the distance of the FS's center from the pivot, R = y / mjj~+n^, 
and on one coordinate, no, while mo is fixed. It is observed that C cr monotonously decreases 
with R, starting from C CI = 2 at R = (we recall again that C& = 2 for the FS in 
the ordinary DNLS equation). Note that there are different pairs (mo, no) which have equal 
values of R = \J m\ + n\ and give slightly different C cr . For instance, for the pair (5,0), 
C CT = 1.51, while C CT = 1.50 for (4,3). Hence, the stability depends on the two-dimensional 
structure of the solution. 

Figures [T] and [3] show that the FS is destabilized, with the increase of C, through the 
appearance of a pair of imaginary eigenfrequencies. Direct simulations of the dynamical 
evolution of unstable FSs in the framework of Eq. (J3j) demonstrate that the instability 
does not destroy the solitary wave. Instead, it transforms the waveform into a persistent 
breathing structure, see a typical example in Fig. [5j 

B. Analytical estimates 

The decrease of C CT with the increase of the distance of the FS from the pivot, R, which 
is the main feature revealed by the above numerical analysis, as shown in Fig. HI can be 
explained using an estimate based on the quasi-continuum approximation. To this end, we 
note that stationary solutions to the underlying continuum equation (OQ) are looked for as 
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FIG. 1: (Color online) Plots of real and imaginary parts (left and right panels) of the eigenfrequen- 
cies of linearization around fundamental solitons in the ordinary non-rotating lattice, with f2 = 
(top panels), and for their counterparts in the present model (with Q = 0.1). In the latter case, 
the soliton is centered at (mo = 3, no = 2). 

= e lAt (j)(x,y), with the function obeying the stationary equation, 

A0 = -C Qv 2 + nLA - e [cos {kx) + cos(%)] $ - $ 3 . (7) 

Here, following the discrete model, we have set £ = v = 0, o = —1, and explicitly introduced 
the spatial-scale parameter, 1/y/C, which is a counterpart of the lattice coupling constant 
in Eq. (J3j). Next, we assume the presence of a soliton with amplitude A and intrinsic size 
I, whose center is located at distance R from the rotation pivot (tantamount to the origin, 
in the present case). First, demanding a balance between terms A0, V 2 0, and 3 in Eq. 

9 




FIG. 2: (Color online) Top panels: Contour plots of the real (left) and imaginary (right) parts of 
a stable fundamental soliton centered at m,Q = 3, uq = 2, for C = 1.5. Bottom panel: the complex 
plane of the stability eigenvalues for this soliton. 

(J7|), and estimating them, respectively, as AA, A/l 2 , and A 3 , we conclude, in the lowest 
approximation, that 

|A| ~ cr 2 ~ A 2 . (8) 

Further, the soliton as a whole will be in equilibrium relative to the rotating lattice poten- 
tial if the action of the centrifugal force, generated by term ~ Q in Eq. (J7]), is compensated 
by the force of pinning to the periodic potential. The estimate of the latter condition yields 
CQR/l ~ e. Substituting here / ~ ^/C/Wl 

as per Eq. (jSJ), we arrive at a final estimate, 
R ~ ^/a/IA]"^ C -1 / 2 , which predicts dependence C cr ~ 1/R 2 (C CT is realized here as the 
largest value of C that can provide for the balance between the centrifugal and pinning forces 
at given R). The latter dependence is qualitatively consistent with the numerical findings 

10 



10 








1 .8 


10 


8 
6 








1 .6 


8 
6 










1.4 




4 










4 










1.2 




2 










2 










1 















c 


-2 








0.8 


-2 


-4 








0.6 


-4 


-6 








0.4 


-6 


-8 








0.2 


-8 


10 










-10 





m 



0.5 
0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 



-0.5 




0.15 



-10 







10 



FIG. 3: (Color online) The same as in Fig. [21 but for an unstable fundamental soliton, placed at 
the same position as in Fig. [21 at C = 1.8. 

showing the decrease of C cr with R in Fig. 0] (at very large C, when the analytical formula 
predicts small R, it is irrelevant, as the above consideration tacitly assumed that the size of 
the soliton was essentially smaller than the distance to the pivot, I -C R', it is irrelevant too 
at very small C, as the quasi-continuum approximation cannot be used in that essentially 
discrete case). 

Using the quasi-continuum approximation, it is also possible to estimate the order of 
magnitude of the rotation frequency, in physical units. In the application to the BEC, 
assuming the lattice spacing ~ 1 /im and the condensate of 7 Li or 85 Rb, which provide for 



31 



32|, 



the possibility of the attraction between atoms and, thus, the formation of solitons 
undoing of rescalings which cast the GPE in the normalized form of Eq. ([1]), we conclude 
that Q = 0.1 corresponds, in physical units, to the rotation frequency ~ 100 Hz or 10 
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FIG. 4: (Color online) Dependence of the stability border for the fundamental solitons, C cr , on the 
distance of the soliton's center from the rotation pivot, R (left panel), and on one of the center's 
coordinates, no, for fixed too (right panel). In the latter figure, the value of uq € [0,8] increases 
progressively from top to bottom. 

Hz, for lithium and rubidium, respectively. As concerns the above-mentioned realization 
of the model in terms of the twisted bundle of optical fibers, an estimate shows that, for 
the carrier wavelength ~ 1 /im and separation between the fibers in the bundle ~ 10 /im, 
Q — 0.1 corresponds to the twist pitch, which we define as a length at which the twist attains 
the angle of 2tt, of the order of 5 cm. 



IV. VORTEX SOLITONS 



Discrete vortex solitons (VSs) of the 2D DNLS equation were systematically developed 



in Ref. 



35| as complex stationary solutions which feature a phase circulation of 2irS around 



the central point, at which the amplitude vanishes, with integer S identified as the vorticity. 
Prior to that, time-periodic multibreather states that may feature a vortical structure were 
found in 2D Hamiltonian lattice dynamical models 36(. The center of the VS may coincide 
with a site of the lattice, or may be located in the middle of a lattice cell; the corresponding 



vortices are called "crosses" (alias rhombuses) and "squares" 



of these states has been studied in detail, both for S = 1 



35 



respectively. T 



le stability 



37(] and S > 1 



29|, E]. In 

particular, the VS crosses with S = 1 (and, as above, with A fixed to be 1) are stable in 
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FIG. 5: (Color online) The time dependence of the squared amplitude of a perturbed unstable 
soliton from Fig. [3] at its center, (mo = 3, no = 2). Notice the robust oscillatory behavior indicating 
the breathing nature of the resulting solution. 

(S=l) 

the corresponding interval of values of the coupling constant, C < C cr = 0.781, and 
the instability above this points transforms the VS into an ordinary fundamental soliton, 
with S = 0. While all VSs with S = 2 are unstable in the ordinary 2D DNLS equation, 
the vortex solutions with S = 3 have their stability interval (more narrow than for S = 1), 



c < a 



(5=3) 



0.198. In all cases, the instability sets in via the Hamiltonian Hopf bifurcation, 



represented by quartets of eigenvalues 

We have constructed localized vortices (of the cross/rhombus type), with S = 1 and 
S = 2, in the rotating-lattice model based on Eq. ([3]), and examined their stability. Unlike 
the above analysis of the FSs, we consider here only on-axis VSs, whose centers coincide with 
the rotation pivot (R = 0). We note that, while the FS with R = is virtually identical to 
its counterpart in the ordinary model, with Q = 0, for VSs the situation is quite different; in 
particular, the VSs centered at the rotation pivot feature quite nontrivial stability properties 
with the increase of rotation frequency Q, see below. Therefore, while the results for the 
FSs were displayed for Q = 0.1 and growing R, in this section we focus on R = but vary 

n. 
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A. Vortex solitons with 5 = 1 



The first noteworthy effect of the rotation on the VS with S = 1 is that, for given Q, its 

( 5=1) 

stability region features not only the upper bound, Ca , as in the usual DNLS model, but 
also a lower one, Cc\~ 1 ^: 

^r i} < c < (9) 

~ ( S=l) 

At a given value of Q, VSs are exponentially unstable for C < C C r , and they feature 
an oscillatory instability, accounted for by a Hopf bifurcation, at C > Ccr ■ This situation 
takes place up at Vt < fiif -1 '* = 0.037. An example is given by Fig. El which displays the 
instability growth rate of the VS, |Im(A)|, versus C for fixed Q = 0.01. The plot also shows 
the real part of the eigenfrequencies, |Re(A)|, and includes, for the sake of comparison, the 
same dependences in the usual DNLS equation, with Q = 0. The figure shows that, in this 
case, Ci 5=1) = 0.31 and C c ( f =1) = 0.788. 

The overall stability region for the VSs with S = 1 in the plane (C, Q) is presented in 
Fig. [71 It is observed that Ccr ~^ slightly increases with f2, while the growth of C& ~^ with 
Q is fast. As a result, at Q > fl C r the stability region does not exist. In the absence of the 
stable VSs, unstable ones feature coexistence of exponential and Hopf instabilities. Figure 
|6] displays an example of the latter regime, and shows that, for Q = 0.04, C K a ' = 0.87 and 
C c ( f =1) = 0.81. 

In fact, the effect of the destabilization of the VSs with S — 1 at small C [C < C cr - , 
as said above], is a higher-order phenomenon, in terms of the expansion of the stability 
eigenvalues in powers of C, when it is assumed small. Indeed, comparing it with the stability 
analysis for the cross vortices in the ordinary model, with Q = 39], we observe the 
following. At Q = 0, modes of small perturbations around the vortex of the cross/rhombus 
type feature a pair of real eigenfrequencies, with A = ±C at the first order in small C 
(in the present notation). The same pair appears in the present context, see Fig. [61 At 
larger C, these eigenfrequencies will give rise to the Hamiltonian-Hopf instability, upon their 
collision with eigenfrequencies bifurcating from the phonon band of linear excitations. As 
shown in Ref. 39J, the DNLS equation with Q = also gives rise to a pair of higher-order 
eigenfrequencies, A ~ ±C 2 in the present notation. The effect of the rotation forces the 
latter eigenfrequencies to separate along the imaginary axis, thus inducing the instability at 
small C . However, for larger C, the pair again becomes real, stabilizing the configuration. 
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FIG. 6: (Color online) Dependences on C of the real and imaginary parts (left and right panels, 
respectively) of eigenfrequencies of small perturbations about vortex solitons with S = 1 in the 
ordinary (non-rotating) DNLS lattice (top panels), and in the rotating one, with 0, = 0.01 (central 
panels) and ft = 0.04 (bottom panels). In the rotating case, the center of the vortex coincides with 
the rotation pivot of the lattice. 
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FIG. 7: The full stability diagram for vortex solitons with S = 1, centered at the rotation pivot. 
The E+H label holds for coexisting Hopf and exponential instabilities. 

Four different regimes identified in the stability diagram in Fig. [7] are illustrated by typical 
examples of the stability and instability of the VSs with S = 1 in Fig. [8] as follows: the 
exponential instability, caused by the higher-order eigenfrequencies, as described above, is 
shown in the top left panel; the top right panel shows a linearly stable case. The oscillatory 
instability is presented in the bottom left panel (the latter case is shown for sufficiently 
large C, to allow the eigenfrequencies, which originally linearly depend on C, as indicated 
above, to collide with eigenfrequencies bifurcating from the phonon band). Finally, the 
mixed oscillatory-exponential instability, which is possible at f2 > f2 cr , is displayed in the 
bottom right panel. 

(S=l) 

Since the destabilization of the VS at C > Ccr is accounted for by the Hopf bifurcation, 
as confirmed by Fig. El this unstable VS is transformed into a persistent breather (not 
shown here), which loses the vortical structure, i.e., is similar to the FS 35[. On the 
other hand, the destabilization at C < Ccr ~ occurs, as seen in Figs. El via a pair of 
imaginary eigenfrequencies with zero real parts, i.e., an exponential instability. Its nonlinear 
development eventually leads to a persistently pulsating localized state with zero vorticity, 
as illustrated by Fig. [91 The transformation to a FS state is also observed in the region of 
the coexistence of exponential and oscillatory instabilities. 
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FIG. 8: Eigenfrequencies of small perturbations around the vortex with 5 = 1 are shown for the 
four different cases: exponential instability for O, = 0.01 and C = 0.2 (top left panel); linear 
stability at O = 0.01, and C = 0.5 (top right panel); oscillatory instability at f2 = 0.01 and C = 0.8 
(bottom left panel); exponential and oscillatory instability at = 0.05 and C = 0.85 (bottom right 
panel). 



B. Vortex solitons with S = 2 



Another way in which the rotating lattice drastically alters the stability features of the 
ordinary 2D DNLS model concerns the VSs with 5 = 2. At Q — 0, all such localized 
vortices are unstable due to an imaginary eigenfrequency proportional (at the leading order) 
to coupling constant C 0, see the top panel of Fig. [10] (which may be compared with 
Fig. 4 of |37|). In the rotating lattice, the solitons with 5 = 2 acquire a finite stability region, 
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FIG. 9: The time dependence of the amplitude of a perturbed unstable vortex soliton with 5 = 1, 
ft = 0.03, and C = 0.4. Note that the point corresponding to this soliton is located in the left 
instability region in Fig. [7J i.e., the vortex is unstable to non-oscillatory perturbations (see also 
the text). 

as manifested by the example displayed in Fig. [TUJ Note that, unlike the stability interval 
for the VSs with 5 = 1, see Eq. Q), only an upper stability border exists for the 5 = 2 
solitons, i.e., the respective stability interval is < C < C a ■ For instance, in the case 

ret 2") 

shown in Fig. [TO, the stability border induced by the rotation is Ccr ~ = 0.12. Note that 
the mechanism of the stabilization of the 5 = 2 VSs in the present model is different from 
that reported in Ref. 



401 ]. where localized vortices with 5 = 2 were stabilized by an impurity 
(inert site) placed at the center. In that case, the unstable eigenmode was suppressed by the 
defect, making all eigenfrequencies real; eventual destabilization occurred due to collisions of 
those real eigenfrequencies with the linear spectrum. Here, the rotation affects the unstable 
(imaginary) eigenfrequency of the 5 = 2 VS, rendering it real for small C . However, as C is 
increased the eigenfrequency eventually becomes imaginary again, leading to the instability 
of the VS. 

The stability diagram for the VSs with 5 = 2 in the (C, fl) plane is presented in Fig. [TT], 
indicating the increasing stabilization effect of larger rotation frequencies. As in the case 
of 5 = 1, the instabilities of the VSs with 5 = 2 transform it into a persistent breather, 



but without the vortical structure. In fact, a similar qualitative cone 
ordinary model, with Q = 0, where all VSs with 5 = 2 are unstable 



usion was made in the 

mm 
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FIG. 10: (color online) The same as in Fig. [6l but for the vortex solitons with S = 2. Top and 
bottom panels correspond to f2 = and Q, = 0.05, respectively. Notice the contrast between the 
absence of instabilities at small C in the bottom right panel and the immediate destabilization in 
the top right panel. 

C. Quadrupole and octupole solitons 

Along with complex solutions for localized vortices with S > 2, the ordinary 2D DNLS 
equation, with Q = 0, gives rise to real solutions in the form of quadrupoles and octupoles, 



that resemble higher-order vortices, but carry no topological charge |25l. |29| . In particular, 
quadrupoles, which include four lattice sites with alternating phases, have their stability 
region, < C < Ccr Uad \ in the model with f2 = 0; the same is true for octupoles, which are 
based on eight sites with alternating phases. In the anti-continuum limit, the quadrupole 



19 




FIG. 11: The stability diagram, similarly to Fig. [7J but for vortex solitons with S = 2. 

and octupole solutions are seeded, respectively, by the following configurations:, 0o,±i — 

l,0±i,o = -1, and 0o,2 = 4>i-i = 4>%\ = 4>-i,o = l,^i,2 = 4>o-i = <j>-i,i = 4>2,o = -1, all 
other sites having </> mj „ = 0. Note that the latter configuration is tantamount to the necklace 
soliton pattern, that was recently observed in a photorefractive crystal with a photoinduced 



lattice [34] 



We have constructed families of quadrupole and octupole soliton solutions in the present 
model with Q > 0, proceeding from the ant i- continuum patterns to finite C. When doing 
so, the rotation pivot in Eq. (jSj) was set at the center of the respective pattern, i.e., we set 
£ = v = and £ = v = 1/2, for the quadrupole and octupole solutions respectively. It was 
found that the critical values of the coupling constant, Ccr"^ and Cc° , which border the 
stability regions for both these families, very weakly depend on Q, unlike the situation for 
the VSs with S = 2, cf. Fig. [Til but quite similar to what was found above for the localized 
vortices with S = 1, see the right stability border in Fig. [7J 

V. CONCLUSION 

The objective of the present work was to introduce a discrete version of the 2D model 
combining the self-attractive cubic nonlinearity and rotating square-lattice potential. The 
discrete model can be implemented in BEC stirred by a rotating strong optical lattice, or, in 
principle, also in a twisted bundle of nonlinear optical fibers. Localized solutions of two types 
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were considered: off-axis FSs (fundamental solitons), with the center placed at distance R 
from the rotation pivot, and on-axis VSs (vortex solitons), with vorticity 5 = 1 and 2. For 
the FSs, the stability interval was found in the form of < C < C CT (R), where C is the 
coupling constant of the discrete lattice, and C CT (R) a monotonically decreasing function, 
see Fig. 0J A qualitative explanation to this result was proposed, based on the analysis of 
the balance between the lattice-pinning and centrifugal forces. For VSs with 5=1, the 
dependence of the stability region on rotation frequency Q was found, in the form of Eq. (Q 
and Fig. [7J with the conclusion that the stability is only possible for Q < Q CT . A key feature, 
which makes the situation different from earlier stability analyses of such structures in the 
standard 2D DNLS model, is a higher-order eigenfrequency, shifted towards instability for 
sufficiently weak lattice coupling, in the presence of rotation. On the other hand, VSs with 
5 = 2, which are always unstable in the model with Q = 0, are stabilized by the rotation 
in region < C < dr~ 2 \ as shown in Fig. [TlJ For the vortices with 5 = 2, the reverse 
effect happens in comparison with 5 = 1, namely an unstable (at Q = 0) eigenfrequency is 
tipped by the rotation in the opposite direction; i.e., for 5 = 1 a real eigenfrequency becomes 
imaginary in the presence of Q, while the reverse is true for 5 = 2. Quadrupole and octupole 
solitons, with the center coinciding with the pivotal point, were briefly considered too, with 
a conclusion that their stability regions are almost the same as in the ordinary model (with 
Q — 0). An estimate for relevant values of Q in physical units was given for both physical 
realizations of the model, i.e., the rotation frequency of the optical lattice stirring the self- 
attractive BEC, and the pitch of the twisted bundle of optical fibers. 

The analysis initiated in this work can be developed in several directions. In particular, 
results obtained in the continuum model considered in Ref. [26( suggest that, at sufficiently 



large Q, the fully localized FSs, with the center shifted off the axis (i.e., continuum coun- 
terparts of the discrete FSs considered in the present work), are unstable (or do not exist), 
while stable ring-shaped solitons may appear instead, with the center of the ring coinciding 
with the pivot. In fact, the shape of these zero- vorticity rings may be similar to that of the 
vortices considered above. 

[181 ] to this considerably more difficult setting. It would be especially relevant to repeat the 
calculation of the eigenfrequencies for the vortices with 5 = 1 and 5 = 2 by means of those 
methods in the presence of rotation, and quantify the impact of Q on the eigenvalues, to 
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rigorously investigate some effects which were outlined above in a qualitative form. 

Another possibility is to consider a discrete limit of the model with the rotating quasi- ID 

potential, such as the one with the repulsive cubic nonlinearity, which was introduced in Ref . 

2t| | . Studies along these directions are currently underway and will be reported elsewhere. 
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